Chaos assisted tunnelling with cold atoms 
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In the context of quantum chaos, both theory and numerical analysis predict large fluctuations 
of the tunnelling transition probabilities when irregular dynamics is present at the classical level. 
We consider here the non-dissipative quantum evolution of cold atoms trapped in a time-dependent 
modulated periodic potential generated by two laser beams. We give some precise guidelines for 
the observation of chaos assisted tunnelling between invariant phase space structures paired by 
time-reversal symmetry. 

PACS : 05.45.Mt Semiclassical chaos (quantum chaos) 
05.60.Gg Quantum transport 

32.80.Qk Coherent control of atomic interactions with photons 
05.45.Pq Numerical simulations of chaotic models 

I. INTRODUCTION 

During the seventies and the eighties, it became gradually clear that classical Hamiltonian chaos profoundly affects 
the temporal evolution and the spectral properties of the corresponding quantum system as compared to the integrable 
case jl]. Some of these features (dynamical localization, scars of periodic orbits etc) share striking similarities 
with concepts originating from condensed matter such as weak and strong localization H. In fact these phenomena 
can be recast in terms of wave transport in disordered media, (quasi)-randomness being of statistical or dynamical 
origin. In this context, it is important to understand the mechanisms underlying a key feature of wave propagation 
which has no classical analog: tunnelling. 

Tunnelling refers to any wave process which is classically forbidden to real solutions of Hamilton equations. For one- 
dimensional (ID) autonomous systems, it is well known that the quantum tunnelling probability through an energetic 
barrier can be evaluated semiclassically with the help of classical complex solutions of Hamilton equations |^|J^| . The 
direct generalization of this procedure to higher dimensional systems is straightforward for separable dynamics but 
is already subtle for still integrable, but no longer separable, dynamics |||7]]. In the generic case of chaotic dynamics, 
it even proves extremely hard to handle and the situation, until recently, seemed hopeless. This is so because, in the 
presence of chaos, the analytical and topological properties of the classical complexified phase space are far from trivial. 
During the last ten years however, theoretical and numerical investigations on autonomous 2D and time-dependent 
ID Hamiltonians systems have started to highlight some mechanisms P-[T2]] and much insight has been gained on the 
influence of such classical non-separable dynamics. Experimental evidence of such mechanisms, which is still lacking, 
would be of great interest especially in the light of the subtle interplay between interferences and disorder. 

In this paper, we consider ID time-dependent dynamics, one of the simplest case where irregular motion can appear, 
and we study chaos assisted tunnelling. Our effective Hamiltonian model, which is derived from an experimentally 
achievable situation, exhibits three main properties. First, its classical dynamics is time-reversal invariant. Second 
it is controlled by a single real external parameter 7 (for 7 = the dynamics is integrable and chaos develops more 
and more in phase space as 7 is increased). Third, there exists in phase space, for a whole continuous range of 7, a 
pair of stable islands X + and X_ which are time reversed images of each other. By stable islands we mean the set 
of regular classical trajectories in phase space which stay near a stable equilibrium point or near a stable periodic 
orbit of the system. In this case, no real classical orbit started in one of these islands can go into the other one. 
However, the quantum dynamics of a wave-packet, initially prepared in one island, will display a periodic behavior. 
The wave-packet jumps from one island to its time reversed image Jl3| . In the quantum spectrum this tunnelling 
process appears via the existence of non degenerate energy doublets whose splitting give the inverse of the tunnelling 
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time between X + and Z_. Varying 7 slowly modifies the geometry of the islands themselves. The crucial point is 
that it will drastically change the classical dynamics for some initial conditions lying between the islands. For 7 small 
enough, the chaotic layers are too small to play a significant role at h scales and hence cannot influence the quantum 
behaviour of the system which is essentially still regular. For larger values, but still before the stable islands are 
completely destroyed, there is a chaotic regime where varying 7 or H (by h we mean here Planck's constant divided by 
some typical classical action) alone induces large fluctuations, on several orders of magnitude, of the doublet splittings 
around their mean value. This in turn means large fluctuations of the tunnelling periods. These wild fluctuations 
induced by small changes of any parameter are a signature of the so-called "chaos assisted tunnelling" regime. It 
has been extensively studied both theoretically and numerically in the situation described above but nas not 

been yet observed in real experiments, the main reason being its extreme sensitivity to small changes in the classical 
dynamics. Any uncontrolled variation of 7, be it noise or dissipation, will dramatically swamp or destroy the signal. 
The observation of this highly fluctuating tunnelling regime thus requires both an accurate control of the dynamics, 
of the preparation of the initial state and of the analysis of the final state. 

Atom cooling techniques provide systems which fulfill all these requirements. They allow an accurate manipula- 
tion and control of internal and external degrees of freedom and are now a useful tool to produce situations where the 
wave character of the atomic motion is essential [ p"5| . A great variety of potentials can be produced to influence the 
atomic motion, be it by means of inhomogeneous magnetic fields, material gratings or laser light. Optical lattices with 
crystalline or quasi-crystalline order flJaTlq] can be easily produced where atoms mimic situations usually encountered 
in condensed matter [ fl9|]2"o(| . Dissipation (spontaneous emission and atom-atom interaction) is easily controlled and 
coherence times of the order of 10 ms can be achieved. This is why cold atoms are a unique tool to study transport 
properties of waves, be it quantum chaos or weak localization p2] , |2^ |. 

The general organization of this paper is the following. In section Owe explain the origin of the effective Hamiltonian 
for the experimental situation under consideration. In section III we study the corresponding classical dynamics and 
show why this effective Hamiltonian is relevant for chaos assisted tunnelling. In section IV we quickly review some 
of the usual theoretical techniques when dealing with both space and time periodic quantum dynamics. We also 
illustrate how some quantum spectral properties have a natural classical interpretation. In section [v] we show, with 
the help of numerical experimen ts, h ow chaos assisted tunnelling arises in our system and then explain how to observe 



it in a real experiment. Section VII is devoted to some concluding remarks. 



II. EFFECTIVE HAMILTONIAN 



A. Light shifts 



The very basic physical mechanism underlying our forthcoming discussion is the following: when an atom is 
exposed to monochromatic light, its energy levels are shifted by the interaction. These level shifts originate from 
the polarization energy of the atom in the incident light field and are called the light shifts |24j]. In the dipolar 
approximation, they only depend on the field intensity value at the center-of-mass position of the atom. If the field 
intensity is space-time dependent, then a moving atom will experience dipolar forces: inhomogeneous light shifts act 
as potentials and alter the center-of-mass motion of the atom. By appropriately tailoring the space-time dependence 
of the light field, one can then produce a great variety of potentials for the external atomic motion as evidenced 
by the atom cooling industry. Note, however, that the atom-light interaction is also responsible for a dissipative 
phenomenon (real absorption of a photon followed by spontaneous emission) which shortens the temporal coherence 
of the atomic wave function. By using a laser light far detuned from any atomic resonance, it is possible to control 
this stray phenomenon and maintain it at a reasonably low rate. 

In the following, we shall describe a simple physical situation for atoms where chaos assisted tunnelling should show 
up. 



B. Experimental configuration 

Although the internal structure (hyperfine Zeeman sub-levels) are of major importance in the atom cooling tech- 
niques, we will here for simplicity model the atom by a two-level system (as indeed, only one optical transition 
usually governs the dynamics). We consider a dilute sample of identical (but independent) two- level atoms propagat- 
ing in the light field configuration created by two monochromatic standing waves with frequencies lo± = lul ± 5w/2 
where 8u <C lu^. We denote by \g) and |e) the ground-state and the excited state of each atom, these levels being 
connected by an electric dipolc transition of angular frequency w a t and width T. All atoms are supposed to be initially 
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prepared in their ground state. Each standing-wave is produced along the x axis by two counter-propagating laser 
beams and we suppose all fields to be linearly polarized along the z axis (see figure 0). After a suitable choice of 
space-time origin, the total electrical field strength is 

E(x, t) = [E + cos(w + t) + E- cos(a>_£)] cos(k L x) (1) 

where E± are the field strengths of the two standing-waves. At this point we have neglected the difference in wave- 
vectors of the standing-waves. For this to hold, it is sufficient to assume that the atomic sample size is small enough. 
Typically, the difference in the k vectors will be of the order of 10 -9 or less (see below), so that this requires the 
atomic cloud to be smaller than typically few kilometers, which is amply satisfied in a standard magneto-optical trap. 

C. Dimensionless effective Hamiltonian 

The effective Hamiltonian which describes the atomic motion is derived in appendix |X| under some common and 
well-controlled approximations. It acts in the Hilbert space of a one-dimensional system which is simply the x 
component (position) of the center of mass of the atom (which means that the internal degree of freedom as well as 
the y and z coordinates can be eliminated, see Appendix). It reads: 

2 

H = - V cos(2k L x)[9 + cos(Sujt)} (2) 
2M 

where Vq = f — M7 + f2_/8c>L and 9 d = (f2+/2f2_) + (fi_/2fi_|_), with 6l = — Wat the detuning with respect to the 
atomic frequency and Sl± — dE±/H (d being the atomic dipole strength). Without loss of generality we will assume Vq 
to be positive since, if Vq is negative, it is sufficient to shift x by n /2fcj, to recover this case. 

def 

In the following, it will prove convenient to work with dimensionless quantities. Rescaling quantities through r = 

Suit, q d = 2kLX, p d = (2kL/M5ui)p x , 7 = f (Ak\/M6u) 2 )VQ and H c g = f (Ak\/M8uj 2 )H , then yields the dimensionless 
effective Hamiltonian: 

P 2 

H cS = — - 7 (6 + cos t) cos q . (3) 

Such an Hamiltonian describes the dynamics of a periodically driven pendulum. The associated quantum canonical 
commutation relation is [q,p] — ih c s and we get h e g = 8ujr/5u: where ujr = %k\/2M is the atomic recoil frequency 
and Slo is the beating frequency between the laser waves. 

Such an effective Hamiltonian clearly exhibits two of the three properties mentioned in the introduction: the 
corresponding classical dynamics is governed by a single classical parameter, the dimensionless coupling strength 7, 
and is invariant under time-reversal symmetry (p, q, r) 1— > (— p, q,— r). It is worth mentioning that the semiclassical 
limit h c g — > is realized here by increasing the beating frequency 5u between the two laser waves. 

With our field configuration, only 6 > 1 can be achieved. As a slight generalization, we extend the range of 9 to 
any positive value since one can design other field configurations where 9 < 1 occurs. For example, 9 = yields the 
Hamiltonian studied in [|l3| in a different context. 

D. Orders of magnitude 

Let us give some typical experimental parameters. For Rubidium atoms, the atomic parameters are M = 85 amu, 
A a t = 2-Kcjijj & t = 0.78 /im, r/27r = 6 MHz, lor/2tt = 3.8kHz and saturation intensity I sa t = 1.6mW/cm 2 . Using 
far-detuned laser beams (28l/T = 10 4 ) focussed down to 500 /im (power lOOmW), with a beating frequency Sui/2tt = 
60 kHz, lead to 7 = 0.4 and heg = 0.05. With such values, spontaneous emission can be neglected up to times of 
the order of few ms. It is worth noticing the tiny energies which come into play (Vq ~ 5neV), by several orders of 
magnitude smaller than the typical ones for mesoscopic systems. 

III. CLASSICAL DYNAMICS 
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A. Poincare surface of section 



A Poincare surface of section provides the usual tool for visualizing the classical dynamics |2q ]. As H c s is 2tt- 
periodic both in time and space, this surface of section simply consists in the whole phase space itself (which has 
the topology of a cylinder) where trajectories (p(r),g r (r)) are seen stroboscopically at every time period 2n. In the 
following, without any substantial loss of generality, we will restrict our analysis to the case 9 — 1 which is easily 
experimentally achieved when the standing waves have the same field strengths. 

Figure H shows stroboscopic plots of phase space orbits for different 7's. For 7 = 0, p is a constant of motion, so that 
the system is integrable and the surface of section is composed of horizontal lines. For 7 weak enough (fig. g-a), the 
orbits remain confined to invariant curves. These invariant curves stratify the whole phase space and the dynamics 
appears regular. One can clearly see well separated stability islands, each being bordered by a separatrix. This is 
the situation encapsulated in the KAM theorem for near-integrable motion: although no globally-defined constants of 
motion exist, some invariant curves can still be constructed which order the dynamics. As 7 is increased (fig. ||-b), 
some of more and more invariant curves are broken and chaotic layers start to spread around separatrices. These layers 
fill some portion of phase space but motion is still predominantly confined to invariant curves. Above some coupling 
threshold (fig. ||-c to ^-e), stochastic orbits invade phase space and the surviving stability islands are surrounded by 
a connected chaotic sea. This occurs for 7 ~ 0.1. The phase space structure in this regime is typical of a mixed 
dynamics where regular orbits co-exist with stochastic ones. If 7 is increased further (fig. |^-f), the stability islands 
disappear (or are too small to be seen at this scale) and one gets global chaos. We note however that, even in this 
situation, the chaotic portion of phase space is still bounded by invariant curves which means that chaos can only 
fully develop within some range of momentum p. 

B. Resonances 

At this stage, let us rewrite the effective Hamiltonian as follows: 



H cS = H + 7-ffi = — 7 cos q - - cos(q -\- t) — cos(q - r) (4) 

The physical interpretation of the various terms is rather simple: the two counter-propagating laser beams at 
frequency uj + create a stationary wave which, in turn, creates for the atom an effective optical potential proportional 
to the square of the modulus of the electric field in the standing wave, hence the cos q dependence (it is actually 
rather 1 + cos q, but the constant term does not play any role in the dynamics). The same effective potential is due 
to the standing wave created by the two lu- counter-propagating beams. A pair of counter-propagating beams at 
frequencies uj + and uj- does not create a standing wave in the lab frame. However, in a frame sliding at constant 
velocity vq = (lu + — w_)c/2w£, (=1 in rescaled units) the two laser beams are shifted in frequency by Doppler effect 
and appear to have equal frequency, building another stationary wave and yet another effective optical potential. In 
the lab frame, this appears as a modulated optical potential moving at velocity vq. By symmetry, there are two such 
effective potentials moving either to the right or to the left. These are the cos(g ± t) terms in the Hamiltonian. 

This form of the Hamiltonian is in order to point out the perturbative terms which may be resonant with the 
unperturbed frequencies. When 7 = 0, the system is integrable since we recover free motion: H c h reduces to 
Hq = p 2 /2 and (p, q) are exact action-angle variables. For 7 > 0, the absence of any constant of motion generates 
chaos. Stroboscopic plots of phase-space trajectories are no more constrained to follow lines of constant Hq but 
generically fill densely a two dimensional volume in phase space. As long as 7 is small enough, these volumes remain 
thin enough not to be distinguished from regular lines at the scale of finite precision of the measurements and/or 
the calculations (c/. figure ||-a). Nevertheless, for higher values of 7, some chaotic layers can be seen (c/. figure y-b) 
between regular regions. They consist of portions of phase space where trajectories are exponentially sensitive on 
initial conditions. From classical first-order perturbation theory |25|, chap. 2], we can infer that a term of the form 
Acos(sq — rr), where (s, r) are integers, will create a resonance of width Ap — A\f~A around the point p — r/s. In our 
case, s = 1 and there exist only three such resonances. They are located at p — (r — 0) and at p — ±1 (r — ±1). 
This can be seen in figure || (a to e). For each resonance there exists one stable and one unstable periodic orbit with 
period approximately 27r|r|. In the stroboscopic plot of the surface of section, they appear as a stable and unstable 
fixed points and give rise locally to the well-known phase-space portrait of a pendulum. In the following, we will 
denote by Xo,X^,X- the three stable islands associated with r = 0, +1, —1 respectively. The physical interpretation 
of these three resonances is simple: each resonance is associated with one of the modulated potentials (either static or 
moving) described above. For example, the fixed point at the center of the 1+ resonance is associated with a periodic 
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orbit where the atom moves at almost constant velocity vq, being is fact trapped in the minimum of the moving 
optical potential. The other two components of the potential appear along this orbit as rapidly varying potentials 
which are adiabatically averaged to constant values. As the atom can be trapped in any of the 3 modulated potentials, 
we obtain three stable periodic orbits at the centers of the 3 resonance islands. 

For 7 small enough, the resonances are well separated and the motion is quasi-integrable. Chaos will develop when 
the resonances start to overlap. This is the celebrated Chirikov's overlap criterion [^6| and its evaluation gives 7 ~ 0.1 
in our case. Thus chaos develops in phase-space regions where the kinetic energy term and the perturbation are of the 
same order of magnitude. Taking into account higher perturbation orders in 7 will shift the position in phase-space of 
the previous resonances as well as the frequency around their stable points. For instance, it can be seen in figure |^-e 
that the stable island X + is centered on a point having a momentum slightly larger than +1. Perturbation terms of 
higher order will also introduce other resonances of smaller size. It is precisely the overlap of the infinite cascade of 
such resonances which gives rise to the chaotic layers. Nevertheless, Chirikov's criterion already gives a good order of 
magnitude for the onset of chaos. For higher 7, the previous three resonant islands of stability have shrunk inside a 
large chaotic sea and eventually disappear completely (c/. figure ||-f). Nevertheless revival of some stable islands can 
still be observed for some narrow windows of high 7's. In our situation, chaos cannot invade the whole phase space 
but is bounded by regular coasts. This is so because chaos develops where resonances overlap. Sufficiently far away 
from the resonances, atoms move so fast that they experience an average time-independent potential. Then chaos is 
absent and one recovers (quasi)-free motion when |p| ^> 1. 

C. Typical classical phase space portrait in the chaos assisted tunnelling regime 

The two resonant islands 2±, when they exist, are related by a discrete symmetry: the time reversal invariance. 
As can be seen in figure ||, the atoms trapped in one island cannot classically escape from it: the boundaries of the 
islands play the role of a dynamical barrier which atoms cannot cross. Hence jumping from one island to the other is 
a classically forbidden process though it is expected to occur in quantum mechanics. This is precisely the tunnelling 
situation we arc interested in. In fact, we will study the tunnelling between X + and X_ for 7 varying from 0.1 to 0.3 
since, in that range, classical chaos may play a revealing role though the two stable islands still occupy a significant 
volume in phase space. 

Note that, in the physical situation described by i? c ff , tunnelling occurs in momentum coordinate instead of space 
coordinate as it is usually presented in standard textbooks. The denomination of "dynamical tunnelling" [^7j, refers 
to this situation. The reason for investigating this situation is that manipulation of cold atoms allows for a better 
control (preparation and detection) of momentum rather than position. 

IV. QUANTUM DYNAMICS 
A. Floquet-Bloch theory 

When an autonomous Hamiltonian is spatially periodic, it is well known |2^] that its spectrum is organized in energy 
bands E n (k). These bands are labelled by a set of integers, the band index n, and depend continuously on a set of 
real numbers, the Bloch numbers k. As E n (k) and the associated eigenfunctions are periodic functions of the k's, all 
the physical information is contained in the first Brillouin zone. For ID-systems, it is simply the interval [— ^, ^[, 
where Q is the spatial period of the Hamiltonian. 

When the Hamiltonian is time-periodic, with period T, the analogous of Bloch theory is Floquet theory |29|~|3l|] . The 
eigenvalues of the evolution operator U(t+T, t) over one period take the form e -leT / ri ° ff . The e's are r-independent real 
quantities which are called the quasi-energies of the system. Due to the time-periodicity, the quasi-energy spectrum 
as well as the associated eigenfunctions are now invariant under e 1— > e + 2irh c ff/T. 

For -ffeffi the application of Bloch and Floquet theorems with Q = T = 2n yields a spectrum made of quasi-energy 
bands e„(fc) where n goes over the whole set of integers (for a detailed derivation see appendix |b|). For brievety, 
we will define \n,k,r) the ket at time r with Bloch angle k, with quasienergy e n (k) and which is a solution of the 

Schrodinger's equation (following the notations of appendix^, we have set |n, k,r) = |V'e„(fe),fc(' r )))- We will also 

define |n, k) = f \n, fc, t = 0). 

As it can be seen in figure [|, the band spectrum has the topology of a torus since it is both periodic in quasi-energy 
(with period h c s) and in Bloch number (with period 1). 
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B. Numerical calculation of the Floquet eigenstates 



As derived in appendix the Floquet eigenstates can be obtained by diagonalization of the Floquet-Bloch operator, 

K 

K(p,q,T, k) = — + ^ cff ^ 7 cos q (1 + cost) - i^eff-r- , (5) 

with periodic boundary conditions both in time and space. The eigenvalues, which depend on the Bloch vector fc, are 
the quasi-energies of the system. The band spectrum is symmetric with respect to the axis k = since the operator (|^) 
is invariant under the transformation k i— » — k and p i— > —p. From the expression of the Floquet-Bloch operator and 
the boundary conditions, it is very natural to expand the eigenstates on a basis set composed of products of the 
type 4>i m (T,q) = exp(inr) exp(imq) which automatically obey the periodic boundary conditions. In such a basis, the 
operator K has very strong selection rules, namely: 

|An| < 1 and |Am| < 1 (6) 

All matrix elements violating one of these selection rules is zero. Hence, the matrix representing the operator K in 
this basis is sparse and banded, and all matrix elements have simple analytical expressions. This is well suited for 
numerical diagonalization (powerful algorithms exist, for example the Lanczos algorithm). All the numerical results 
presented here use this method. We checked that the effect of the truncation of the basis is negligible: the size of 
the Floquet matrix is considered as sufficiently large when increasing it modifies the value of the quasi-energy on the 
scale of the numerical noise only, say 10 -15 in double precision. Not only this criterion is a proof of algorithmical 
convergence but also it is a safeguard against numerical discrepancies since we are looking for exponential small 
quantities. 



C. Husimi representation and classification of the quantum states 

Classical dynamics is very illuminating when describing the states \n,k,r). In order to strengthen the corre- 
spondence between classical phase-space structures and quantum states, it is convenient to work with the Husimi 
representation of quantum states [p2| . 

Such a representation associates to each quantum state \ip) a phase space function ?/> H (p, q) (where p and q are real 
numbers) defined by 

r(p,q)=N^ \(zM 2 (7) 

where \z) is the normalized coherent state corresponding to the complex number z = (q + ip)/\/2h c g. is a (p, q)- 
independent normalization factor. Because \z) is a minimal gaussian wave packet with average momentum p and 
average position q, the Husimi function ip H (p, q) contains some information about the degree of localization of \ip) in 
phase space. 

The minimal cell size in phase-space allowed by the Heisenberg inequalities is ft e ff- Let us see how classical phase 
space structures of typical size larger than h e f{ are mirrored at the quantum level. In figure |5| and ^ we have plotted 
some of the e n (k) corresponding to the Hamiltonian (||) for specified fixed values of 7 and h c g. Von Neumann and 
Wigner arguments |33| claim that, generically, no exact degeneracy can occur: one rather gets avoided crossings. Of 
course, this is relevant provided the minimal energy splitting is greater than the resolution in energy. Some Husimi 
functions are plotted in figure || and || (a to /) []. 



1 A technicality should be mentioned here: |u e ,fc) and \ip t ,k) obey some spatial boundary conditions which are lost when working 
with their Husimi representations, essentially because the coherent states do not fulfill themselves these properties. To deal 
with spatially (quasi-)periodic phase-space functions it is necessary to unfold the coherent states ^] into 

\z,k) d = J2 eimkQ \ z + m Q) ( 8 ) 
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In appendix it is shown that we have 

v n .k = 7f / (», fc, t] p \n, k, r)d,T = " , (10) 

which generalizes the velocity theorem p8| , Appendix E] to periodic time-dependent hamiltonians. 

Figure ||-b shows an example of the Husimi representation of a state with sufficiently high average velocity to be alike 
a free eigenstate of Hq. Far from quasi-degeneracies it is localized in a narrow strip of width Ap ~ 2irh c s/ Aq <~ h c ff 
(since Aq covers 2n) and which is centered on one of the two classical phase space trajectories of energy about k /2 
(compare with figure ^-e). Its quasi-energy band (figure |^) is an arc of the parabola of the free motion but can hardly 
be distinguished from a straight line of slope v„,fc if k is restricted to one Brillouin zone. We will naturally call these 
states quasi-free states. 

Some states have their Husimi functions localized in the resonant stable islands (in Xq but also in T± ) . The number of 
these states is semiclassically given by the volume of these islands divided by 2nh c s. Far away from quasi-degeneracies, 
these states are at any time centered on the stable periodic orbit: it can be explained within a semiclassical approach 
and can be observed in figure ^-c for a state localized in lo and in figures [§-d,e for states in Z_, X + respectively. Their 
average velocity as well as their Husimi functions depend exponentially weakly on the Bloch parameter k. 

The last class of states which can be encountered corresponds to chaotic ones, that is to say states whose Husimi 
functions are negligible on a typical distance of \Jh c s out of the chaotic seas (c/. figure ||-a). Unlike the previous 
ones, their Husimi functions are very sensitive to any variation of k since they are delocalised in the whole chaotic 
sea which spreads over all elementary cells. The large classical distribution of possible velocities is to be linked to the 
very fluctuating slopes of the quasi-energies as k is varied. 

D. Tunnelling states 

Although the system as a whole is of course time-reversal invariant, this is no longer true for its restriction at 
a fixed value k of the Bloch vector. Indeed, the operator K is not time-reversal invariant, because of the crossed 
term kp. In other words, the time-reversed partner of a state with Bloch vector fc is a state with Bloch vector — fc. 
It is only at the special value fc = (and also fc = 1/2 since fc is defined modulo 1) that K is invariant under time 
reversal symmetry. Therefore it corresponds to the typical situation of tunnelling between X + and Z_ . Every island 
state localised about \p\ = 1 is quasi-degenerate with another one. These doublets represent a symmetric and an 
antisymmetric combination of states localised in one island only (c/. figure ^-f). The energy splitting Ae n of these 
states for fc = is precisely the signature of tunnelling: it is nh c ff divided by the typical time an atom takes to jump 
from one island to its time-reversed image, i.e. to reverse the sign of its velocity. 

V. CHAOS ASSISTED TUNNELLING 
A. Large fluctuations 

After having selected the two quasi-energy bands corresponding to the two states which are localized the more 
deeply inside the islands T±, we are able to plot the splitting as a function of h c g . The great advantage of studying 
fluctuations when the effective Planck constant is varied is that it does not affect the classical dynamics. The behavior 
of the splitting is very different whether chaos is present at the fi e ff scale or not. In the chaotic regime (c/. figure |?] 
and ||), that is when h c g varies in a range where chaotic seas can be resolved, the splittings vary rapidly versus 
the change of any parameter, in our case h Q g. Moreover, the variations of the splittings, despite being perfectly 
deterministic, are apparently erratic - without any regular structure - and cover several orders of magnitude. They 
show that direct coupling to the chaotic sea is the key mechanism for their understanding and are a signature of chaos 
assisted tunnelling pa,M. In fact these huge fluctuations are reminiscent of the universal conductance fluctuations 



and define for instance 

< fc (p,g,r) d = N\(z,k\i, n , k (r))\ 2 (9) 
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observed in mesoscopic systems [|37],[38| since tunnelling is nothing else than wave transport from one stability island 
to the other. In contrast, in the regular regime where chaotic seas are smaller than h e s, the splittings are expected 
to vary smoothly jf| . 

In figure |?], we show the splittings of the pair of states localized at the center of the resonances islands I± , as a 
function of h c s. They display huge fluctuations over about 4 orders of magnitude while the general trend is a decrease 
as h e s — ► 0. Similarly, when plotted as a function of 7, see figure ||, they also display large fluctuations. The general 
trend is here a fast increase of the typical splitting with 7; this is associated with the shrinking of the regular island 
when 7 increases which results in a increasingly large tunnelling probability The overlaps of the regular states (still 
supported by the islands) with the chaotic states increases. Therefore the coupling between the two components of 
the tunnelling doublets which involves the chaotic states increases as well. 

In order to understand both the general trend and the origin of the fluctuations, two points of view can be used: a 
quantum point of view and a semiclassical one. 

B. Semiclassical interpretation 

If the chaotic sea is large, it is rather intuitive that it can be easier to first tunnel from the center of the regular 
island to the chaotic sea, then propagate freely in the chaotic sea to the vicinity of the symmetric island and finally 
tunnel to the center of the symmetric island than directly tunnelling between the two islands. The crucial point is 
that because the chaotic component is explored rapidly and densely, it does not cost anything to cross the chaotic 
sea. Tunnelling trajectories can be viewed as complex trajectories (i.e. with complex position and momentum) which 
are real at both the starting (in the initial island) and ending (in the symmetric island) points. The tunnelling 
amplitude associated with a single tunnelling orbit is essentially cxp(— lm(S) / h c g) where S is the complex action of 
the tunnelling orbit. In a usual 1-dimensional system (like a double- well) , there is only one such trajectory at each 
energy and the tunnelling rate thus displays the well-known exponential decrease. In a chaotic system, it may happen 
that there is a whole set of tunnelling trajectories whose imaginary parts of the action are essentially identical. In 
such conditions, the actual tunnelling amplitude is the sum of all individual amplitudes (each taken with its proper 
phase) which results in a very complicated quantity which fluctuates when parameters are changed. In some sense, 
this is analogous to the speckle pattern obtained when plenty of optical rays with various geometries are randomly 
interfering. This is the very origin of the (deterministic) fluctuations of the tunnelling rates and consequently of the 
energy splittings. The general trend (exponential decrease) is related to the typical imaginary part of the action of 
the tunnelling trajectories. 

C. Quantum point of view 

A complementary quantum point of view is possible. One can divide the eigenstates of the system in two subsets: 
the "regular" states localized in the resonance islands and the "chaotic" states localized in the chaotic sea. The two 
sets are only weakly coupled by tunnelling. Because there are two symmetric islands, the regular states are essentially 
doubly degenerate (neglecting direct tunnelling). The chaotic sea also has the two- fold symmetry and states can be 
classified as even or odd. The two series of odd and even states are ignoring each other. Hence, when by accident, a 
even chaotic state is almost degenerate with a even regular state, they repell each other: at the same time, there is 
usually no odd chaotic state with the same energy. Thus, the odd regular state is not significantly repelled. Hence, 
the splitting appears because of different shifts of the even and odd regular states. Close to any avoided crossing 
between either the odd or the even regular state and a corresponding chaotic state, a large splitting is obtained. On 
the contrary, far from any avoided crossing, the splitting is small. Hence, the fluctuations are just associated with 
the existence of a large number of successive avoided crossings. The typical size of these fluctuations is related to 
the typical size of the avoided crossings while the typical parameter range of these fluctuations is the distance (in 
parameter space) between two consecutive avoided crossings. A model implementing this idea (each regular state 
is independently and randomly coupled to the chaotic states of the same symmetry) has been proposed in |j9[ and 
further used in (4C[ |. In this model, the chaotic states are modeled by a Hamiltonian belonging to the Gaussian 
Orthogonal Ensemble of random matrices while the coupling between the regular state and the chaotic state is also 
taken as a random Gaussian variable. With these assumptions, the splitting distribution can be calculated. Let us 
denote by A the mean level spacing between chaotic states and by a the typical strength of the coupling between 
the regular states and the chaotic sea. Only if a <C A is the regular state weakly coupled to the continuum (if this 
inequality is violated, the regular state is completely diluted in the chaotic sea by the strength of the coupling) . We 
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thus assume the inequality to be valid. Then, the distribution of splittings Ae is given by: 



where 



P(Ae) = - „ \ , for|Ae|<a . . 

v ' ir s 2 + Ae 2 (11) 
P(Ae) ~ for |Ae| > a 



s = ^~ (12) 



The interpretation is rather simple. The maximum splitting is observed exactly at the avoided crossing where the 
levels are shifted by ±cr/2 on both sides of their unpertubed positions. Hence, the splitting cannot be larger than 
about a. In fact, there is an exponentially decreasing tail in the distribution P(Ae) (associated with the Gaussian 
fluctuations of a) which we do not detail here because it is not relevant in our present case, s is the typical splitting 
one expects to observe: it corresponds to the shift typically due to the closest chaotic state. The full distribution is a 
(truncated) Cauchy distribution: it is obtained as the overall effect of all chaotic states lying above or below in energy. 
Note that, in the absence of the truncation of the Cauchy distribution, the average splitting is not defined because the 
corresponding integral diverges. Hence, it is better to discuss the typical splitting s rather than the average splitting. 
It is the slow decrease of the Cauchy distribution for large splitting which is responsible for the huge fluctuations in 
the splittings, which can be as large as a ^S> s. This is reminiscent of random processes such as Levy flights where 
rare events are dominant 

In fig. |p, we show the statistical distribution of splittings that we obtain numerically in the chaotic regime (normal- 
ized to the typical splitting in order to work with s = 1). The distribution is shown on a double logarithmic scale and 
compared with the Cauchy distribution. On can clearly see two regimes: for small |Ae|, P(|Ae|) is almost constant 
and decreases with a slope -2 for large |Ae|. The agreement with the Cauchy distribution is very good, which proves 
that the model catches the essential part of the physics in this system. 

The typical splitting s is proportional to the square of the tunnelling matrix clement from the initial state to 
the chaotic sea. Hence, it is expected to decrease roughly like exp(— 2Im(S)//j e g) where S is the complex action of 
tunnelling orbits (the mean level spacing scales as a power of ft e ff and is thus a correction to the main exponential 
decrease). This is roughly what is observed in fig. |7[ Note however that there is a plateau in the range 15 < h~g < 30. 
A similar observation has been done in jl2|, fig. 3] . Although this is not a crucial problem (the statistics of the splitting 
distribution is not affected), a detailed explanation of this behaviour is still lacking, see however p2[ . Finally, it should 
be interesting to calculate explicitely some of the complex tunnelling orbits in our specific system in order to compare 
the imaginary part of their actions with the slope in fig. ^. This work is currently under progress. 



VI. EXPERIMENTAL SIGNAL 

As it can be seen in figure the splittings we want to measure correspond to very tiny scales among the other 
structures in the band spectrum. For say h e g ~ 0.1, the tunnelling times are about h e g/Ae ~ 10 3 times the typical 
period (T~ few/xs). The observation of atoms having reversed their velocity is therefore still possible since during 
1 ms spontaneous emission has not begun to spoil our dynamical model. Nevertheless, measuring tiny splittings which 
are hidden so deeply in the spectrum is far from being straightforward. Several steps are needed: preparation of the 
initial state, the experiment itself (where chaos assisted tunnelling takes place) and the analysis of the final state. 
During the first step, one must prepare a state localized inside one resonance island, i.e. localized both in position 
and momentum spaces. The idea is to use an adiabatic transfer from an initial state extended in position space by 
slowly branching the effective potential. Although the second step looks trivial (one just has to wait), the dispersion 
in the vecto r ang le k makes things much more difficult as only the k = states are related by time-reversal symmetry 



(see section IV D ) and thus tunnel relatively fast. It is however possible to overcome this difficulty as explained below. 
Finally, the detection should be rather easy, using velocity dependent Raman transitions. We now explain in detail 
how the various steps can be worked out. 



A. Adiabatic preparation of the atoms in one lateral stable island 

The first step consists in preparing an initial cloud of cold Rubidium atoms in order to have it located in one of the 
stable island, say T+ , only. Using a standard magneto-optical trap, one can obtain a more or less thermal distribution 
of atoms with velocity of the order of few times the recoil velocity w roc — HI^l/M = 6 mm/s. However - as shown below 



9 



- this is probably too much for a g ood measurement of the tunnelling splitting. Additional techniques (side-band 
cooling |fl3f , Raman cooling [[4~4| , |l5| ) make it possible to obtain a sub-recoil velocity distribution, i.e. atoms with an 
average momentum po and a thermal dispersion Ap x — M Av — Mav Tec with a significantly smaller than 1 . We chose 
the initial momentum to be M5uj/21tL so that, on average, the atoms exactly follow one of the sliding standing wave 
created by a pair to± of laser beams. 

The next step is to slowly (i.e. adiabatically) switch on the standing waves. During this phase, the spatial periodicity 
is preserved, and the Bloch vector k is thus a conserved quantity. Initially, the momentum p x is nothing but the Bloch 
vector (modulo a integer multiple of the recoil momentum). Thus, by preparing a sub-recoil initial state, one populates 
only a small range of k values and, for each k value populated, a single state (momentum eigenstate). Otherwise 
stated, the initial momentum distribution has become a statistical mixture of Bloch states with Afc = a/2 in a single 
energy band. More generally, if the initial momentum distribution is not a sub-recoil one, but has a width equal to 
a recoil momenta (with a > 1), about a bands will be populated. 

Switching on the standing waves increases the optical potential Vq and therefore 7 from zero and enlarges the 
resonance islands. In the sliding frame moving with velocity Po/M = Suj/2kL, the atoms feel a pendulum like 
potential — ^7COsg (in scaled coordinates) in addition to some rapidly time varying terms. Consequently, they will 
adiabatically localize in the potential minima, that is at the center of the resonance island. Increasing 7 localizes 
successively an increasing number of states in X + . The switching time must be sufficiently long as compared to the 
beating period (in this case the discarded terms are still rapidly oscillating) but also to the inverse of the minimum 
energy gap (of the order of h~g if Afc is sufficiently narrow). In order to trap all the initially populated states, 7 has 
to be sufficiently large. For a given Bloch angle, we want to localize the a first states. The quantum energies of a 
pendulum are given by the eigenvalues of the Mathieu equation jl6|, chap. 20] (see figures [l0| and 11 ) . For a pendulum 
whose hamiltonian is 

2 

ffpcnd = — - - cos q , (13) 

the phase space volume enclosed by the separatrix is given 161/7/2. Semiclassically, it corresponds to 16^^/2/ (2irh c g) 
states. This number will be of the order of a when 7 reaches the value: 

(aTrficff) 2 

7adiab ^ ■ ^ > 

Figures [l(] and [ll] show the exact energy levels of the system together with the ones using the pendulum approxi- 
mation and the Mathieu equation, as well as plots of selected Husimi representations for few eigenstates. This allows 
to check that : 

• The pendulum approximation works well in the regime of interest, say up to 7 = 0.1. 

• The semiclassical estimate of the number of trapped states is sufficiently accurate for our purpose. 

• The Husimi representations of the trapped states are well localized: especially, the states in figure [ll] have two 
well separated components in the 2+ and I_ islands, meaning that we are in a real case of tunnelling. 

• During the initial increase of 7, the "ground" state is well isolated (in energy) from the other ones, which means 
that an adiabatic preparation is possible. 

Nevertheless, for this adiabatic preparation to be valid, we must make sure that 7 has not reached a range where 
chaos have non negligible effects on quantum properties. Physically, we want that chaotic layers at 7 = 7adiab have 



a small volume compared to the Planck constant. From figure 10, we can obtain the upper bound limit of 7 by 
estimating when the (small) avoided crossings become too large to be passed diabatically. For 1/50 < h c a < 1/5, we 
observe that chaos have no influence for 7 < 7 c haos — 0.04. 

The adiabatic preparation of all the atoms in one stable island will be achieved if 7 a diab < 7chaos that is if the atoms 
are cold enough to have 



a < ~ — . (15) 

7Tft e ff tleS 

In every situation considered in the following, we will have to check that this condition is fulfilled. Next we have 
to reach the desired value for 7 (~ 0.18) in the chaotic regime while preserving the state in the island. This will 
be achieved if 7 is increased sufficiently fast so that all encountered avoided crossings with chaotic states are passed 
diabatically. 

Hence the whole preparation of the initial atomic state proceeds by two steps: a first adiabatic increase of 7 at the 
very beginning to significantly populate regular states in one island followed by a diabatic increase of 7 to preserve 
them in the chaotic regime. 
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B. How to force tunnelling ? 



As discussed in section IV D| , the tunnelling splitting is small only for Bloch vector k = 0. As it is not presently 



possible to prepare only this value of k in a real experiment, it seems at first sight that only a small fraction of atoms 
(close to k = 0) may effectively tunnel, hence considerably reducing the signal to noise ratio. A solution is to force 
all atoms to go through the k ~ region. The simplest idea is to impose a slow increase of the k value, by adding a 
constant external force F from outside. Then, the full Hamiltonian which governs the dynamics is: 

H'(p,q,T) = H(p,q,T)-qF . (16) 

The potential V induces a dynamical drift in the Bloch angle: 



(17> 

It is shown in appendix |c| that this relation, which is well-known in the time independent case (see for instance |47| , 
chap. 6]), remains valid when H is periodic in time. 

A convenient way of realizing experimentally such a constant force is to chirp the laser frequencies, that is make all 
the frequencies drift linearly in time. In an accelerated frame, the laser frequencies appear as constant, and we are 
back to our model. However, in this non incrtial frame, the constant acceleration is translated in a constant force, 
hence the system is governed by equation (|l6|). This method has been used with cold atoms, see |^9| . 

The global result is a slow drift of the k distribution. This makes that the various k classes successively come close 
to k — and are thus able to tunnel. Whether the atom will effectively tunnel or not depends on the time scale on 
which k changes. If k varies rapidly, the avoided crossing at k = is crossed diabatically, i.e. the velocity distribution 
will not be modified. If k varies slowly, it is crossed adiabatically (rapid adiabatic passage). The Landau-Zener 
formula Eal y ields for the typical time scale for the crossover between diabatic and adiabatic crossing. 

Figures|l2| and [l3| illustrate, for theg — 0.2037 the drift of atoms initially localized in 2+ (one energy band) whose 
distribution in k covers 1/10 of the Brillouin zone from 1/10 to fc ~ (sub-recoil initial velocity distribution 

with width of the order of v iec /5i) Let us have the atoms evolving under the force F in order to get a global translation 
of 1/10 in k. After a sliding across the splitting at k = 0, if the force F is weak to follow the energy level adiabatically, 
the average momentum of the atoms has reversed its sign as can be seen on figure Measuring the critical value 
of F for this Landau-Zener like transition to occur furnishes a mean to measure the splitting. 

Another possibility would be to modulate the external force (and thus the k values) periodically in time in order 
to induce a resonant transfer between X + and I_ . 

In any case, the method may work only if no other avoided crossing come into play. Numerical investigations show 
that there are mainly tiny avoided crossings along the energy curve of interest. However, as a general rule, there 
are usually few avoided crossings with similar or larger size than the avoided crossing of interest. If such an avoided 
crossing is also passed adiabatically, it will of course spoil the momentum distribution. Hence, it is crucial for the 
initial k distribution to be sufficiently narrow to avoid this problem. Hence, a sub-recoil velocity distribution seems 
necessary. 



C. Detection 



After the atoms have interacted with the modulated waves, one can switch off the lasers either abruptly or adi- 
abatically (in which case the atoms adiabatically leave the resonance island). In both cases, the atoms which have 
tunneled from X + to Z_ will end with a momentum close to — po- A standard time of flight technique should be enough 
to detect them. A more sophisticated technique, for example based on velocity sensitive Raman transitions [^9 43 
could also be used with a subrecoil resolution if needed [|19| . 



VII. CONCLUSION 



In this paper we have proposed a simple and accessible experimental configuration in which the observation of 
chaos assisted tunnelling should be feasible. It consists in atoms propagating in the light field of two far-detuned 
monochromatic standing waves with slightly different frequencies. Observing the tunnelling effect however requires 
sub-recoil cooling techniques to conveniently prepare the atomic sample together with a well-controlled experimental 
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procedure (adiabatic preparation of the atomic state in one stable island followed by a drift of the Bloch vector) . As 
the tunnelling period fluctuates over several decades when the potential strength varies (see fig. 8), observing these 
fluctuations requires a stabilization of the laser intensity at the level of a few percent. 
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APPENDIX A: DERIVATION OF THE EFFECTIVE HAMILTONIAN 

To derive the effective Hamiltonian (|J), we basically proceed in three steps |5(| : 

- We first assume that any dissipation process can be safely ignored. Indeed for tunnelling to be observable, 
phase coherence of the atomic wave-function must be preserved during all the process. In our case, this means 
that spontaneous emission must be negligible. It can be shown |5l| ] that increasing the laser detuning considerably 
decreases phase coherence loss. Hence the evolution will be essentially Hamiltonian provided each laser beam is 
sufficiently far-detuned from the atomic resonance : 

|<J±| = |w±-w at |»r (Ai) 

When this holds, the total Hamiltonian operator for this two-level system is the sum of three terms: the kinetic 
energy operator describing the center-of-mass motion of atoms of mass M the energy operator for the internal degrees 
of freedom and the coupling between internal and external degrees of freedom. In the dipolar approximation, the 
interaction is just dE(x, t) and does not depend on y and z. The dynamics along y and z is just trivially described 
by free motion and can be easily eliminated since the total quantum state factorizes as a plane wave in y and z. One 
is then left with the dynamics along x which is described by the hamiltonian : 

flat = ||(l e X e l + \9M) + ^at|e)(e| - dE(x,t)(\e)(g\ + \g)(e\) (A2) 

where p x is the atomic momentum along x. \g) and |e) are the ground-state and excited-state and d is the atomic 
dipole strength connecting them. 

- Second, we expand the total atomic state as: 

|*) = ip g {x,t)\g) +tp e {x,t)exp(-iuj L t)\e) (A3) 

and we drop high frequency (optical) anti-resonant terms as far the amplitudes change slowly during an optical period. 
This is known as the rotating wave approximation [52| and features the averaging procedure to eliminate fast variables 
in classical perturbation theory |25|. This yields the coupled amplitude equations: 



We = -^l^e ~ hS L i> e - ^L^A^ g . (A4b) 

In these equations Sl = — w a t is the mean laser detuning, the star denotes complex conjugation and Q(x,t) 
reads: 

Q(x, t) — [f2 + exp(— iduj t/2) + exp (idut/2)] cos(k L x) (A5) 

where Cl± = dE±/h are the Rabi frequencies of each standing wave. 

- As a final step, we assume now that atoms initially prepared in their ground-state mostly evolve in their ground- 
state. This means that the whole atomic dynamics is solely determined by the ground-state amplitude ip g . For this to 
hold, adiabatic eli mina tion of the excited-state amplitude must be justified. If the spatial partial derivatives were 



absent, equations (A4) would just describe the Rabi oscillation phenomenon. It is then known that far off resonance, 
i.e. when the frequency separation of the states is much larger than any other frequencies, the Rabi oscillation is very 
small in amplitude. A sufficient condition is : 
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\6 L \>n ± ,6u> (A6) 

If in addition we assume that the excited-state kinetic energy is very small (which will be easily achieved with cold 
atoms): 



\5l\ » i>, 



^ (A7) 



2M 

then adiabatic elimination of the excited-state amplitude just amount to neglect the spatial and temporal derivatives 



of ip e in equation (A4t) which is then solved as ip e ~ — {Q,*/28i)ipg <C ip g . It is then easy to see that the ground-state 
amplitude ip g obeys an effective Schrodinger equation with Hamiltonian: 

2M 4S L K ' 

Eventually, up to an irrelevant purely time-dependent term, we get 

V 2 

H = f2- - V cos(2k L x) [9 + cos(fci)] (A9) 
where V d = -W.+Q-/8SL and 6 d = (fi+/2fi_) + (fi_/2fi+). 

APPENDIX B: FLOQUET-BLOCH FORMALISM 

In this appendix, we briefly recall the Floquet-Bloch formalism which is used for a quantum problem whose Hamil- 
tonian H is periodic both in space and in time. We will denote T and Q the temporal and spatial periods respectively. 
Let us first consider the time periodicity. We define 29 3l|]53|| 



K(p 7 q,T) d ^ ~ih-^-+H(p,q,T) (Bl) 

where p and q stand for canonical Hcrmitian operators whose commutator is [p, q] = —ih. 

If U(t',t) denotes the unitary evolution operator from r to r' associated with Hamiltonian H, the periodicity of 
the dynamics implies that U(t + T,T) = U(t, 0) and 

U(t + T,t) = U(t + T,T) [/(T,0) U(0,t) = [t/(0, 7-)]" 1 U(T,Q) U(0,t) . (B2) 

This shows that U(t + T, t) and U (T, 0) differ by a unitary transformation and hence have the same spectrum (but 
of course different eigenvectors), independent of r. The eigenvalues of U(t + T,t) have unit modulus and can be 
written as e~ ltnT / h where e n is the so-called quasi-energy defined modulo 2-kH/T. If |-0„(r)) denotes the corresponding 
eigenvector, we can define the Floquet state: 



\Xn{T))=e lenTin \Mr)) (B3) 

which is by construction periodic with period T . 

Inserting the definition of the Floquet state in the time-dependent Schrodinger equation, we immediately obtain: 

K(p,q,r)\xn(r)) = e n \xn(r)) . (B4) 

which means that the quasi-energy spectrum is obtained by diagonalizing the Floquet Hamiltonian in the space of 
time-periodic functions. 

The second step consists in making use of the invariance of H under spatial translations with period Q. The unitary 

translation operator Tq = e _1 P^/ R commutes with K. We can then use the spatial counterpart of the Floquet theorem, 
namely the Bloch theorem p8|] , and label the eigenstates of K with the Bloch number k 6 [— n/Q, n/Q[ (the first 
Brillouin zone), which means diagonalising K in each subspace with fixed k. If one defines the Floquet-Bloch states 
as: 



13 



KfcM) - e-^| X „, fc (r)) = e^^/h-^Mr)) , (B5) 

where {ipn,k{T~))} forms a complete orthogonal eigenbasis, it is easy to show that they can be obtained by diagonalizing 
the Floquet-Bloch Hamiltonian: 

K(p,q,T,k)=K(p + hk,q,T) (B6) 

on the subspace of time and space periodic functions. In our specific case, the Floquet-Bloch Hamiltonian reads: 

( v + hk) 2 d 
K(p, q, r, k) = — — — 7 (d + cos r) cos q - ih— . (B7) 

2 GET 

The spatial periodicity of the Floquet-Bloch states leads to a discrete set of dispersion relations e n (k) . For fixed n, 
the set of all quasi-energies e n (k) for k in the first Brillouin zone [—tt/Q,tt/Q[ is called the nth band of the system. 

Let us now obtain the velocity theorem (|To|). Using the above relations, we have (ip et k(' T )\p \'*Pe,k{ T )) = {u ei k(r)\ (p + 
Hk) \u e ^( T )} — ( u e,k( T )\(h~ 1 dK /dk^\u e k(r)) . The derivation with respect to k of the relation (u ej fe(r)| K |w £) fe(r)) = 
e(k) leads to 



(^e,fc(r)|p|^ e ,fc(r)) 

1 d 



hdk \Ue,k( T )\ ) ^|Me,fe(r)) 



(B8) 



The two last terms of the the right hand side are opposite since the normalization of the u's leads to d(u e ^k\u €j k) /dk = 0. 
Moreover, after time-averaging (B8) over T, the total r-derivative vanishes since the u's are precisely T-periodic while 
the time-independent fc-derivative of the quasi-energy remains unchanged. Eventually, 

^j\^Ar)\pHeAr))dr = ^. (B9) 

APPENDIX C: BLOCH ANGLE DYNAMICS 

In this appendix we derive equation ( |l7| ) which is valid for an arbitrary strength of the constant force F provided 
that the potential V = —Fq remains strictly linear in q. Let us choose a state \ip(r)) evolving under H 1 = H + V 
such that it coincides with a Floquet-Bloch state at t = 0: 

iheS ^)l = H'(p,q,T)\ip(T)) (CI) 

and 

|V(t = 0)) = \i(j n , k (r = 0)) . (C2) 
In the interaction picture we have immediately 

ifteff d| ^ (r)) = -FUHt, 0)qU(r, 0) IV' (t)) (C3) 

CLT 

where = f W {t,Q)\^i{t)) and where U denotes the evolution operator under H. Let |0(r)) be the ket defined 
by 

\<t>( T )) d = e- iF ^/^«|^(r)) . (C4) 

It is straightforward to obtain its evolution : 

d^r)) = G(r) mr)) (C5) 
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where 




Since W(t,Q) commutes with the translation operator Tq, it can be checked that [G(t),Tq] = 0. The evolution 
of \4>ij)) under G will therefore preserve its initial quantum number k : 




which shows that |V>( r )) is actually a Bloch wave with a Bloch angle given by fc(r) = k(0) + Fr/hes even if it spreads 
among the quasi-energy bands. 
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FIG. 1. Experimental configuration under consideration: a cloud of two- level atoms is exposed to two monochromatic 
standing waves with frequencies lo± — lol ±8lu/2 (Slo -C ^z,). All fields are linearly polarised along the same direction and are 
sufficiently far-detuned from the atomic resonance so that dissipation effects can be ignored. 
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FIG. 2. Stroboscopic plots of trajectories inphase space for different initial conditions at r = and different 7's. The 
classical dynamics is governed by Hamiltonian (pi) with = 1. At low 7 values, resonance islands are visible separated by 
quasi-free motion. As 7 increases, the resonance islands grow and chaos appears close to the separatrices. The situation of 
interest for chaos assisted tunnelling is when two symmetric islands are separated by a chaotic sea, as 1+ and X_ in (d) and 
(e). 
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FIG. 6. Quasi-energy bands for 7 = 0.18 and ties = 0.2037 and some Husimi representations of typical states. It is a 
zoom of the dashed squared zone in figure U For k / 0, one can find states, like in d or e, whose Husimi function is localised 
in one stable island I_ or 2T+ (compare with figure fye). I^the frame where the center of the island is fixed, these states 
correspond to the "ground" states (some excited states may of course exist if h e g is small enough as can be seen in figures 
[12] and [l3|). For k — 0, we recover time reversal symmetry through the existence of quasidegenerate doublets of symmetric or 
antisymmetric combinations. This is the typical tunnelling situation: following adiabatically with k the state in e (which has 
an average velocity about +1), by decreasing k we get the state in d which has a reversed velocity. This reversal of the velocity 
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FIG. 7. Fluctuations of the energy splittings Ae„ between pairs of symmetric/antisymmetric states localized in the X± 
resonance islands (here shown for the "ground state" inside the island). The classical dynamics is fixed at 7 = 0.18 (compare 
with figure ^-e). The existence of large fluctuations over several orders of magnitude is a signature of chaos assisted tunnelling. 
On the average, In |A e „ | appears to decrease more or less linearily with except for the plateau at 15 < < 30, 
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FIG. 8. Fluctuations of the energy splittings Ae n between the pair of symmetric/antisymmetric states (H e g is fixed at 19.309) 
as a function of 7. Again, the large fluctuations over several orders of magnitude are a signature of chaos assisted tunnelling. 
The global increase with 7 is due to the growth of the chaotic sea as 7 increases, see figure 0. 
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FIG. 9. Statistical distribution of the energy splittings Ae (normalized to the typical splitting) between pairs of symmet- 
ric/antisymmetric states localized inside the X± resonance islands (7 = 0.18, k = 0), represented on a double logarithmic scale. 
On can clearly distinguish two regimes: constant at small Ae followed by a I/Ae 2 decrease and finally a rapid cut-off (not 
shown in the figure). The solid line is the Cauchy distribution predicted by Random Matrix Theory. 
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